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Abstract: This paper studies two stochastic bulk demand (S, s) inventory models A and B both with phases 
PH demand occurrence and k2 phases PH lead time distributions respectively. In the models the maximum 
storing capacity of the inventory is S units and the order for filling up the inventory is placed when the inventory 
level falls to s or below. The demands occur when the absorption occurs in the demand arrival PH process. Its 
size is random and is bounded above with distribution function depending on the phase from which the 
absorption occurs. Lead time for an order realization has phase type distribution with the supply size > S and is 
bounded above with bounds depending on the phase from which the absorption occurs in the PH lead time 
distribution. When the inventory has no stock, demands wait for supply forming a queue. Only after clearing the 
waiting demands, the inventory is filled up. When an order is placed only after its realization, the next order can 
be placed. When the inventory is filled up due to an order realization, units received in excess are returned. In 
model A, the maximum of the maximum demand sizes in the phases is greater than the maximum of 
thek2supply sizes. In model B the maximum of the maximum demand sizes in the kj^phases is less than the 
maximum of the k2supply sizes. Matrix partitioning method is used to study the inventory systems. The 
stationary probabilities of the inventory stock size, the probability of the waiting demand length, its expected 
values, its variances and probabilities of empty levels are derived for the two models using the iterated rate 
matrix. Numerical examples are presented for illustration. 

Keywords: Block Partition, Bulk Demand, Matrix Geometric Method, Phase Type Distribution, Supply after 
Lead Time. 



I. INTRODUCTION 

In this paper two bulk demand (S, s) inventory systems are studied with phase type (PH) distributions. 
Models with PH distributions are very useful and important since PH distribution includes Exponential, Hyper 
exponential, Erlang and Coxian distributions as special cases and serves as suitable approximations for arbitrary 
general distributions as noted in Salah, Rachid, Abdelhakim and Hamid [1]. Inventory and queue models have 
been analyzed by many researchers. Thangaraj and Ramanarayanan [2] have studied two ordering level and unit 
demand inventory systems using integral equations. Jacob and Ramanarayanan [3] have treated (S, s) inventory 
systems with server vacations. Bini, Latouche and Meini [4] have studied numerical methods for Markov 
chains. Chakravarthy and Neuts [5] have discussed in depth a multi-server waiting model. Gaver, Jacobs and 
Latouche [6] have treated birth and death models with random environment. Latouche and Ramaswami [7] have 
studied Matrix Analytic methods in stochastic modeling. For matrix geometric methods and models one may 
refer Neuts [8]. Rama Ganesan, Ramshankar and Ramanarayanan [9] have analyzed M/M/1 bulk arrivals and 
bulk service queues under varying environment. Fatigue failure models using Matrix geometric methods have 
been treated by Sundar Viswanathan [10]. The models considered in this paper are general compared to existing 
inventory models. Here at each demand epoch, random numbers of units are demanded and the maximum 
number of units demanded may be different in different phases. When there is no stock in the inventory, after 
the lead time, realized orders can clear various number of waiting demands. Usually bulk arrival models have 
partitions based on M/G/1 upper-Heisenberg block matrix structm-e with zeros below the first sub diagonal. The 
decomposition of a Toeplitz sub matrix of the infinitesimal generator is required to find the stationary 
probability vector. Matrix geometric structures have not been noted as mentioned by William J. Stewart [11] 
and even in such models the recurrence relation method to find the stationary probabilities is stopped at certain 
level in most general cases indicating the limitations of such approach. Rama Ganesan and Ramanarayanan [12] 
have presented a special case where a generating function has been noticed in such a situation. But in this paper 
the partitioning of the matrix with blocks of size, which is the maximum of the maximum number of demands in 
all phases and the maximum of the order supply sizes in all phases together with PH phases, exhibits the matrix 
geometric structure for the (S, s) bulk demand inventory system with PH distributions. The (S, s) inventory 
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systems of M/M/1 types with bulk demands, bulk supply and random environment have been treated by Sundar 
Viswanathan, Rama Ganesan, Ramshankar and Ramanarayanan in [13]. In this paper two models (A) and (B) of 
bulk demand (S, s) inventory systems with PH demand, PH lead time and infinite storage spaces for demands 
are studied using the block partitioning method to obtain matrix geometric results. In the models considered 
here, the demand sizes are bounded discrete random variables with distinct distributions corresponding to the 
phase from which the PH arrival process moves to absorption state. The lead time distribution is of PH type and 
the order is realized at the epoch at which the absorbing state is reached. The size of the supply is finite 
depending on the last phase before the absorption. When the inventory becomes full, units realized in excess are 
returned immediately. Always waiting demands are given priority and to be cleared before providing stocks for 
the inventory. When the waiting queue of demands is longer than the supply size then the entire supply is 
utilized for reducing the queue length. Model (A) presents the case when M, the maximum of all the maximum 
demand sizes is bigger than N, the maximum of order supply sizes. In Model (B), its dual case, N is bigger than 
M, is treated. In general in waiting line models, the state space of the system has the first co-ordinate indicating 
the number of customers in the system but here the demands in the system are grouped and considered as 
members of M sized blocks of demands (when M >N) or N sized blocks of demands (when N > M) for finding 
the rate matrix. The matrices appearing as the basic system generators in these two models due to block 
partitions are seen as block circulants. The stationary probability of the number of demands waiting for service, 
the expectation, the variance and the probability of various levels of the inventory are derived for these models. 
Numerical cases are presented to illustrate their applications. The paper is organized in the following manner. In 
section II the (S, s) inventory system with bulk demand and order clearance after the lead time is studied with 
PH distributions in which maximum M is greater than maximum N. Various performance measures are 
obtained. Section IE treats the situation in which the maximum M is less than the maximum N. In section IV 
numerical cases are treated. 

II. MODEL (A) MAXIMUM DEMAND SIZE M > MAXIMUM SUPPLY SIZE N 

2.1 Assumptions 

(i) The time between consecutive epochs of bulk arrivals of demands has phase type distribution {a_, T) where T 
is a matrix of order with absorbing rate Tq = —Te to the absorbing state k^+l. When the absorption occurs, 
the next bulk demand arrival time starts instantaneously from a starting state as per the starting vector a = 

(ai,a2., ttfcj) and HfiiCtj = 1- Let (p be the invariant probability vector of the generator matrix (T + TgO). 
When the absorption occurs in the PH arrival process due to transition from a state i to state fci +1, Xi number 
of demands arrive with probabilities P (Xi= j) = Pj for 1 < j < M; and X^iPj=lwhere Mj forl< i <kiis the 
maximum size. 

(ii) The maximum capacity of the inventory to store units is S. Whenever the inventory level falls to s or below, 
orders are placed for the supply of units for the inventory. Arriving demands are served till the inventory level 
falls to 0 after which the demands form a queue and wait for order realization. 

(iii) The lead time distribution of an order has phase type distribution (P , S) where S is a matrix of order k2 

with absorbing rate Sq = —Se to the absorbing state fc2+l and the staring vector is /? = (P^, /?2- , Pi^^) where 
Xfi]^ Pi = Let (|) be the invariant probability vector of the generator matrix (5 + S^P). During the lead time of 
an order, another order cannot be placed. When absorption occurs due to a transition from state i to state ^2+1, 
an order is realized with the supply of constant N, > S units for 1< i < fe2 and the waiting demands are served 
first and the balance if any becomes stock for the inventory. In case the inventory is filled up, the units which 
are in excess, if any are returned immediately. After the realization of an order if the inventory level becomes s 
or below s or the inventory is still empty with or without waiting demands, then the next order is immediately 
placed. When n demands are waiting for 0 < n < N-S at the order realization epoch, n waiting demands are 
cleared, the inventory is filled up and units in excess are immediately returned. When the waiting number of 
demands n at the order realization epoch is such that N-S < n < Nj all the n demands are cleared and N-n units 
become stocks for the inventory if N-n < S and if N-n > S the inventory is filled up and the units in excess are 
returned. If n > Nj demands are waiting when an order is realized, Nj demands are cleared reducing the waiting 
demand length to n-Nj. 

(iv) The maximum of the maximum demand arrival size M=maX]^ <; M; is greater than the maximum of the 
order realization size N=maX]^ <, Nj . 

2.3Analysis 

The state of the system of the continuous time Markov chain X (t) under consideration is presented as follows. 
X(t)={(k, i) : for 0<k< S-s-1 and 1 < i < k^)}\J {(0, k, i, j) ; for S-s<k<M-l; l<i<k^-l <j <fe2}U{(n, k, 
i,j): for 0<k<M-l; 1 < i < fei; 1 < j < ^2 and n > 1 }. (1) 
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The chain is in the state (k, i) when the number of stocks in the inventory is S -k for 0 < k < S-s-1 and the 
arrival phase is i for 1 < i < fej^. The chain is in the state (0, k, i, j) when the number of stocks in the inventory is 
S-k for S-s < k < S without any waiting demand, the arrival phase is i for 1 < i < fe^ and the order-supply phase 
is j for 1 < j < fe2. The chain is in the state (0, k, i, j) when the number of waiting demands is k-S for S+1 < k < 
M-1, arrival phase is i for 1 < i < fe^ and the order-supply phase is j for 1 < j < fe2 The chain is in the state (n, k, 
i, j) when the number of demands in the queue is n M H- k- S, for 0 < k < M- 1 and 1 < n < co, arrival phase is i for 
1 < i < fei and the order-supply phase is j for 1 < j < . When the number of demands waiting in the system is r 
> 1, then r is identified with the first two co-ordinates (n, k) where n is the quotient and k is the remainder for 
the division of r -H S by M; r = M n -H k - S for r>l,0<n<oo and 0 < k < M- 1 . Let the survivor probabilities of 
arrivals be P( > m ) = = 1 - I]„™ iPn, for 1 < m < M; -1 with Pq' = 1- for all i, 1 < i < (2) 
The chain X (t) describing model has the infinitesimal generator of infinite order which can be presented in 
block partitioned form given below. 

-fii So 0 0 . . . • 
B2 Ao 0 . . . • 



Qa= 



0 

Aa 



0 



0 

Aa 



0 



(3) 



In (3) the states of the matrix are listed lexicographically as 0, 1, 2, 3, .... For the partition purpose the states in 
the first two sets of (1) are combined. The vector 0 is of type 1 x [k^^S — 5) -I- fe]^fe2(M-SH-s)] and n is of type 
1 x(/Cife2M). They are 0 = ((0,1),(0,2). ..(0,fei)... (S-s-1, 1)... (S-s-1, /<:i),(0,S-s,l,l),(0,S-s,l,2)...(0,S-s, /Ci,/C2) 
...(0,S,1,1)...(0,S, fei,fe2),(0,SH-l,l,l)...(0,S-i-l,/ci,fc2)...(0,M-l,l,l)...(0,M-l,feifc2)). For n > 0 the vector is 
n=((n,0,l,l),(n,0,l,2)...(n,0,l,fc2),(n,0,2,l),(n,0,2,2)...(n,0,2,fc2),(n,0,3,l)...(n,0,fei,fe2),(n,l,l,l)...(n,l,/Ci,fe2),(n 

,2,1,1) (n,2,/Ci,/<:2)....(n,M-l,l,l), (n,M-l,l,2) ...(n,M-l,l,/<:2),(n,M-l,2,l) (n,M-l,/<:i, /<:2)). The 

matrices Biand A^ have negative diagonal elements, they are of orders fei(5 — 5) -I- fe]^/c2(M-S-i-s) and kik2M 
respectively and their off diagonal elements are non- negative. The matrices Aq andA2 have nonnegative 
elements and are of orders fe]^fe2M. The matrices Bq and B2 have non-negative elements and are of types 
[ki(S - s) + feife2(M-SH-s)] X (fei/c2M) and (feifc2M) x [k^^S - s) + feifc2(M-S-i-s)] respectively and they are 
all given below. Let © and® denote the Kronecker sum and Kronecker products respectively. 
Let 2i=TeS = {mi,^) + (I,^®S) (4) 
where I indicates the identity matrices of orders given in the suffixes and is of order feife2. Let Tq = 
(tg, tg, ... tg^)' be the column vector of absorption rates in PH arrival process. Let Sg = (Sq,Sq, ...Sg^)' be the 



0 

column vector of absorption rates concerning the PH lead time distribution. Let 

Toj = {tlpj, tlvj tl'v^^y fori < i < M (5) 

where toPj is the rate of absorption from state i to state k-^+l when j demands occur for 1< i < fe^ and for 1 < j < 
M. Let the matrix Aj = [Tgj a ] for 1 <j <M, (6) 

S,j = {s'oSj_,^,si Sj_,^ fori < i < N ; S < N; < N and 1 < i < k2 (7) 

where 5; j = 1 if i = j and 5; j = 0 if i ^ j, 1 < i, j < N. Sgj is a 0 column vector for j 7^: iV; for 1 < i < ^2 . Sq^j. is 
a column vector with only one non-zero element (5'o«j)j=5o for 1 < i < fc2 and (5'o«j)y=0 if i T^j for 1 < i, j < fe2- 
Let Uj = 4j (g) [SgjP^ for l<j <N with order k^k2; 
U'j = 4^ ®[Soj]for l<j <N with type k^k2X k^; 
A'j =]^Qj a ®p ] for l<j<M with type k^x (/C1/C2) ; 
Vj = Ylj+.U'j, forl<j<N-l. 



A"j =[Toj a] fori <j <S-s 



1 with order k^, 
(8) 
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(13) 



In B2, the first fei(5' — 5) + /C]^fe2(M-S+s-N) are 0 columns. In (14) the structure presents the case when M-N- 
S+s > 0. When M-N < S-s, then in the column blocks from M-N to S-s all Uj may be replaced by U'j . 
Similarly in the matrix the column blocks from 2 to S-s all Uj may be replaced by U'j. The basic generator 
of the bulk queue which is concerned with only the demand and supply is a matrix of order k^k^M- given below 
in (17) where 2^ =^0+^1+^2 (15) 

Its probability vector w gives, wg^i =0 and w. e = 1 (16) 
It is well known that a square matrix in which each row (after the first) has the elements of the previous row 
shifted cyclically one place right, is called a circulant matrix. It is very interesting to note that the matrix 2^ = 
+ Ai+ A2 is a block circulant matrix where each block matrix is rotated one block to the right relative to the 
preceding. 
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In (17), the first block-row of type kik2X Mkik2 is, W = (Q^ + A^,A-^, A2 , ...,/1m_«_2, -^m-w-i- -^^M-iv + 

-^^-2 + f^2' -^M-i + ^1) which gives as the sum of the blocks (2i ^i+-'^2~^---~^^m-n-2 ~^ 
Am-n-1 + ^M-N + U^+---+Am-2 + U2+ ^M-i + Ui= (T+Too) e(S+5o^) = (T+roa)(8)42 + IkM^ + 
5o^) whose stationary vector is (p®(p. This gives (p®(p{Qi + -4^) + <p(8)(/>I]?l"i'^~^-4i+ i??®*^ 'Tli=i{-^M-i + 

= 0. So {<p0(p, (p<S>(p ■■; <p0(p ). W = 0 = {<p0(p, (p(S>(p ■ ■ ^(plSxfi) W. Since all blocks, in any block-row 
are seen somewhere in each and every column block structure (the matrix is block circulant), the above equation 

shows the left eigen vector of the matrix 2a is {(p®<p,(p0(f), , (p®^)). Using (16), this gives probability 

vector f^^^ (jg^ 

\ M M M M J ' 

Neuts [8], gives the stability condition as, '•n Aq e < w A2 e where w is given by (18). Taking the sum cross 
diagonally in the Aq and A2 matrices, it can be seen that the stability condition can be simplified as follows. 

=^ Y,n=i n (p [ro„a_ ]e (g) (pe =^ i:n=i n Sj'ii (Pj t^oPn= (Pj t'oE(j(j) < w ^2 e = ^ (p<S>(t)Qln=i nUJe = 

^ Q:L,ni(p®(t>Uje = ^ a:Lini(p®(t>Xlk,®SonP)e = ^ ZLin cp e ® cpS^nPe 

^ Hn=i ^ '^'jti 0; ^0 ^n,Nj = ^ Syli 4>j ^0 ^ where ipiand (pj are components of cp and ^ respectively for 1 < i < 

ki and 1 < j < fe2- So the inequality for the steady state reduces to Xfii <Pi to^iXi) < Hfii 4>i ^o^i (19) 
This is the stability condition for the bulk demand (S,s) inventory system with PH distributions when the 
maximum demand size in all arrival phases is greater than the maximum supply size in all supply phases. When 
(19) is satisfied, the stationary distribution of the queue length of waiting demands for units exists Neuts[8] 
Let 31 (k, i) for 0 < k < S-s-1 and 1 < i < fci; ti (0, k, i, j) for S- s < k < M-1, 1 < i < fc^and 1< j < ^2; n (n, k, i, j), 
for 0 < k < M-1, 1 < i < fei, 1 < j < fe2 and n > 1 be the stationary probability of the states in (1). Let tTq = (7i(0, 

1) ...7r(0, fci)...rt(S-s-l, l)...7i(S-s-l, fci) 7i(0, S-s, 1, 1), 7i(0, S-s, 1, 2)...7i(0, M-1, /Ci,/C2)) be of type 
1 x [fci(S -s)+ kj^k2(M-S+s)]. Let 7r„= (7i(n, 0, 1,1), 7i(n, 0,1, 2) ... 7i(n, 0, l,fci), 7t(n, 0, 2, 1), 7i(n,0, 2, 

2) ,...3i(n,0,fci,fe2 ), 7i(n,l,l,l), 7i(n, M-1, 1,1), 7i(n, M-1,1,2). . .7i(n, M-l,feife2)) be of type 1 x/(:ife2Mforn> 

1. The stationary probability vector n = {nQ,ni,n2, ■■■) satisfies the equations 7rQ^=0 and 7ie=l. (20) 
From (20), it can be seen that n^Bi + niB2=0. (21) 
noBo+niA^+n2A2 = 0 (22) 

" 0: for n > 2. (23) 
Introducing the rate matrix R as the minimal non-negative solution of the non-linear matrix equation 

^o+R^i+^?^^2=0, (24) 

it can be proved (Neuts [8]) that 7r„ satisfies the following. 7r„ = for n > 2. (25) 

Using (21), TTg satisfies tTq = 71^52 (26) 

So using (22) , (26) and (25) the vector can be calculated up to multiplicative constant since satisfies the 

equation tt^ [B2 (-Si)"^Bo + + ^^2] =0- (27) 

Using (20), (25) and (26) it can be seen that [B2 i-B^y'^e+(l-R')~'^e] = 1. (28) 
Replacing the first column of the matrix multiplier of tt^ in equation (27), by the column vector multiplier of 
TTi in (28), a matrix which is invertible may be obtained. The first row of the inverse of that same matrix is Ui 
and this gives along with (26) and (25) all the stationary probabilities of the system. The matrix R is iterated 
starting with/?(0) = 0; and by finding R(n + l')=-AoAj^-R'^(n)A2Aj^, n > 0. (29) 
The iteration may be terminated to get a solution of R at a norm level where ||i?(n -I- 1) — 7?(n)|| < £. 
2.3. Performance Measures 
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for r > 0, P (Inventory level is r) = P (INV= r) =] 



(1) The probability of the demand length L = r > 0, P (L = r), can be seen as follows. Let n > 0 and k for 
0 < k < M-1 be non-negative integers such that r = nM + k-S. Then using (21) (22) and (23) it is noted that 
P (L=r) = Xfii ii n (n, k, i,j ), where r = n M + k -S > 0. 

(2) P (waiting demand length L = 0) = P (L = 0) = YtX^ Sfii n (k, i) -nXfii ^i=s-s (0, k, i,j ) and 

Sfii 7r( 5 - r, 0 for s + 1 < r < S. 

(Xfii Y-li TT (0, 5 - r, t,; ) , for 1 < r < s 
P (Inventory level=0, demand length L =0) = Xfii Sj'ii (0,5, i,j ) 

(3) The expected demand length E (L) can be calculated as follows. Demand length L = 0 when there is stock in 
the inventory or when the inventory becomes empty without a waiting demand. So for E(L) it can be seen that 

E(L)=i:t"/+i Xfil Zjl, 7r(0,k,j,i)(k-S)+i:^=i Sfo' Sfii I^jli n in, k, (nM+k-S). 

Let 5i= (0,0,---0,l,l,...l,2,2,...2,...,M-l-S,M-l-S,...,M-l-S)'be a column of type [(S-s) fei+(M-S+s) k^k2\ x 1 
and in the vector, the number 0 appears [(S-s) ki+{?>+\) kik2] times, and the numbers 1,2,3,.., (M-l-S) appear 
kik2 times one by one in order. The vector (0,0,. . .0,1,1,. .. 1,2,2,. . .2,. . .,M-1, M-1, M-1)' where the all 
numbers 0 to M-1 appear feife2 times. On simplification using equations (25)-(28) 
E (L) = TTo^i +M7ri(I-R)-2e + 7ri(I-R)-i52- S7ri(I-R)-ie (30) 

(4) Variance of the demand length can be seen using VAR (L) = E (L^) - E(L)^. Let ^3 be column vector 
53 = [0,..0,l2,...l2 2^..2^... (M- 1-5)^... (M - 1 -5)2]' of type [(S-s) fei-i-(M-SH-s) /<:ife2] x 1 where the 
number 0 appears [(S-s) fe]^-i-(s-i-l) /C]^fe2 ]times, and the square of numbers 1,2,3,.., (M-l-S) appear fei/c2 times 
one by one in order and let ^4 = [0,..0,1^... 1^ 2'^,..2'^, ... (M - 1)^ ... (M - 1)^]' of type Mk^k2xl where 
the number 0 appears feifc2 times, and the square of the numbers 1,2,3,.., (M-1) appear kik2 times one by one 
in order. It can be seen that the second moment, E(L^)=Y,'i^S^^i Xfii ?J'(0,k,i,j)(k-S)2 + 
Jln=iTik=o T!^iiT!^=iT[ (n,k, [M n + k —S]^ . Using Binomial expansion in the second series it may be 
noted E(L2)=7ro53 -I- n(n - l)7r„ e -I- 'Zn=i nn^e]+ 'Zn=i T^n + 2M Y.n=i n n^S2 
-2S Y.n=i I.k=o Ifii (n, k, {Mn + k) Y.n=i Y.f=a Sfii ^ ^' ^'j) ■ After shnphfication, 
E(L2)= nQSj,+M^[n^{l - Ry^2R e + n^il - Ry^e]+ni{I - Ry^S^ + 2M n^il - Ry^S2 

-2S [M7ri(I-R)-2e-i-7ri(I-R)-i52] +S^n^(J - RY^e (31) 
Using (30) and (31) variance of L can be written. 

(5) The above partition method may also be used to study the case of (S,s) inventory system in which the 
supply for an order is a finite valued discrete random variable by suitably redefining the matrices for 1 < j < N 
as presented in Rama Ganesan, Ramshankar and Ramanarayanan for M/M/1 bulk queues [9] and for PH/PH/1 
bulk queues [14]. 

III. MODEL (B) MAXIMUM DEMAND SIZE M < MAXIMUM SUPPLY SIZE N 

The dual case of Model (A), namely the case, M < N is treated here. (When M =N both models are applicable 
and one can use any one of them.) The assumption (iv) of Model (A) is changed and all its other assumptions 
are retained. 



3.1. Assumption 

(iv). The maximum demand size M=maX]^ <; Mi is less than the maximum supply size N=maX]^ Nj . 

3.2. Analysis 

Since this model is dual, the analysis is similar to that of Model (A). The differences are noted below. The state 
space of the chain is as follows presented in a similar way. The state of the system of the Markov chain X(t) is 
X(t) = {(k, i) : for 0 < k < S-s-1 and 1 < i < fci)} U {(0, k, i, j) ; for S-s < k < N-1; 1 < i < fei; 1 < j < /C2} U{(n, k, 
i, j): for 0 <k <N-1; 1 < i < fci; 1 < j < /<:2 and n > 1}. (32) 
The chain is in the state (k, i) when the number of stocks in the inventory is S -k for 0 < k < S-s-1 and the arrival 
phase is i for 1 < i < fej^. The chain is in the state (0, k, i, j) when the number of stocks in the inventory is S-k for 
S-s < k < S without any waiting demand, the arrival phase is i for 1 < i < fej^ and the order-supply phase is j for 
1 < j < /C2. The chain is in the state (0, k, i, j) when the number of waiting demands is k-S for S+1 < k < N-1, 
arrival phase is i for 1 < i < fc^ and the order-supply phase is j for 1 < j < fe2.The chain is in the state (n, k, i, j) 
when the number of demands in the queue is n N H- k- S, for 0 < k < N-1 and 1 < n < 00, arrival phase is i for 1 < 
i<ki and the order-supply phase is j for 1 < j < fe2. When the number of demands waiting in the system is r > 1, 
then r is identified with the first two co-ordinates (n, k) where n is the quotient and k is the remainder for the 
division of r -H S by N; r = N n H- k - S for r > 1, 0 < n < 00 and 0 < k < N-1. The infinitesimal generator Qg of the 
model has the same block partitioned structure given for Model(A) but the inner matrices are of different types 
and orders with distinct inner settings. 
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In (33) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... For the partition purpose the states 
in the first two sets of (1) are combined. The vector 0 is of type 1 x [k^^S — s) + fe]^fe2(N-S+s)] and n is of type 

1 x(/Cife2N). 

0=((0,l),(0,2)...(0,/ci)...(S-s-l,l)...(S-s-l,fei),(0,S-s,l,l),(0,S-s,l,2)...(0,S-s, fei,fe2)...(0,S,l,l)---(0,S, fei./ca), 
(0,8+1,1,1) ... (0,S+l,fei,fe2)...(0,N-l,l,l)...(0,N-l,feife2) ) andforn>l, n= ((n,0,l,l), (n,0,l, 2) ... (n, 
0,l,fe2),(n,0,2,l),(n,0,2,2)...(n,0,2,fc2),(n,0,3,l)...(n,0,fei,fe2),(n,l,l,l)...(n,l,fei,fe2),(n,2,l,l)...(n,2,fei,fc2).... 
(n,N-l,l,l),(n,N-l,l,2) . . .(n,N-l,l,/c2),(n,N-l,2,l). . .(n,N-l,fci, /(:2)). The matrices B\andA\ have negative 
diagonal elements. They are of orders — s) + /C]^/c2(N-S+s) and fe]^fc2N respectively and their off diagonal 
elements are non- negative. The matrices A'q andA'2 have nonnegative elements and are of order fe]^fe2N. The 
matrices B'q and B'2 have non-negative elements and are of types [k^^S — s) + fe]^fe2(N-S+s)] x (/C]^fe2N) and 
(fe]^fc2N) X [ki{S — s) + fc]^fe2(N-S+s)] respectively and they are given below. Using Model (A) for definitions 
ofyly A'j andA"j,for 1 <j <M,andUj, U',,Vj fori <j <N and letting 2i=T©S, the partitioning 

matrices are defined as follows. 
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(34) 
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The matrices B'gand B'2are similar to i4'Q and i4'2. The zero matrices appearing in the first S-s row blocks of 
B'q are of type fci x ( fcife2) and after the S-s th row block the matrices blocks are of order kik2- The model 
presented in (37) and (38) considers the case when M < N-S+s and the blocks A\ falls inside B^.When M > N- 
S+s, then from the row block i+1, /I'y falls outside B\ and A'j is to be placed in B'oup to the row block S-s in 
place ofyly where i=S-s-NH-M , without changing other terms in (34) .The matrix B'2 is similar tOv4'2.In the 
first S-s column blocks U'j appears instead of Uj . 
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(39) 



The basic generator which is concerned with only the demand and supply is Qg = A'^ + A\ + A' 2 and is 
presented in (40). This is also block circulant. Using similar arguments given for Model (A) it can be seen that 

11-1- ■ f<p®<P (P®<p (p®(p {p®(p\ , , 1 ... . , 

Its probability vector is \ ' ™ stability condition remains the same. Following the 

arguments given for Model (A), one can find the stationary probability vector for Model (B) also in matrix 
geometric form. All performance measures including expectation of demands waiting for supply and its variance 
for Model (B) have the form as in Model (A) except M is replaced by N 
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(40) 



IV. NUMERICAL CASES 

For numerical illustration it is considered that the demand arrival time PH distribution has representation 

^ ^ r-2 7 1 

with a = (.3, .3, .4) and the service time PH distribution has representation S= ^ 



T= 



7-4 7 
2 7 -5J 
with p = (.4, .6). 

Six examples are studied. The maximum demand size and maximum supply size are fixed as M=8 and N=8 in 
examples 1 and 2. Examples 3 and 4 treat the cases M=8, N=7 and examples 5 and 6 treat the cases M=5, N=8. 
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The order of the rate matrix R is 48 since it is the product kik2M and kik2N depending on M >N or N > M. The 

probabilities of bulk demand sizes and bulk supply sizes are presented as follows in the examples. 

In the examples 1 and 3 the bulk demand probabilities of sizes 1, 5 and 8 are {pi i=.5, Pi^=.3, Pi g=.2), 

(P2,i~-4' P2,5 =-4, P2,8=-2), (P3,i=-6, P3 5=.3, P3 g =.1) in demand phasesl, 2 and 3 respectively. 

In the examples, 2 and 4 the bulk demand probabilities of sizes 1, 5 and 8 are {pi i=.5, Pi^=A, pis=.l), 

(P2,i~-'^' P2,5 =-5, P2,8=-l)' (P3,i~-6' P3,5=-4i P3,8 =0) in demand phasesl, 2 and 3 respectively. 

In example 5, the bulk demand probabilities of sizes 1 and 5 are {pi i=.5, Pi^=.5, Pi g = 0), (P2,i=-4, P2,5=-6, 

P2,8 — 0), (P3 i=.6, P3 5 =4, P2 8 = 0) in demand phasesl, 2 and 3 respectively. 

In the example 6, the bulk demand probabilities of sizes 1, 5 and 8 are {pi i=.6, Pi^=A, Pi 8=0), (P2,i=-5, 
P2,5=-5, P2,8=0)' (P3,i=-7' P3,5=-3, P3,8= 0) in demand phasesl, 2 and 3 respectively. In all the examples, Pij = 
P2 = P3 =0 for j = 2, 3, 4, 6 and 7. 

The bulk supply sizes in examples 1, 2, 5 and 6 are iV; = 8, for PH phases i=land 2 respectively. The bulk 
supply sizes in examples 3 and 4 are iV; = 7, for PH phases i= 1 and 2 respectively. 

The iteration for the rate matrix R is performed for the same 15 number of times in all the six examples and the 
performance measm^es are written using the matrix iterated R(15) matrix of order 48. The difference norms of 
convergence are presented in table 1 along with expected demand lengths and the variances obtained for the 
examples. The inventory stock level probabilities, probability of both stock level and demand queue length are 0 
and various block level probabilities are presented for the examples 1 to 6 in the table 1 . 

The performance measures obtained show significant variations depending on M, N and pi j . The two levels 
namely the inventory levels and block demand queue levels are both important in the inventory models. Their 
probability values are presented in figures 1 and 2. 



Tablel.Results obtained for six examples 
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Figure 1 Probability of Inventory Levels 



Figure! Probability of Block Levels 



V. CONCLUSION 

Two (S, s) inventory systems with bulk demand and bulk supply after the lead time have been treated. 
Identifying the maximum of the demand size and supply size for forming the demand blocks of such maximum 
sizes along with PH phases, the stationary probability vectors are presented when the inter arrival times bulk 
demands and lead time distributions for supply have PH distributions. Matrix geometric solutions have been 
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obtained by partitioning the infinitesimal generator by grouping of demands and PH phases together. The basic 
system generators of the bulk systems are block circulant matrices which are explicitly presenting the stability 
condition in standard forms. Numerical results for (S,s) inventory systems by varying bulk sizes are presented 
and discussed. Effects of variation of rates on expected queue length and on probabilities of queue lengths are 
exhibited. The PH distribution includes Exponential, Erlang, Hyper Exponential, and Coxian distributions as 
special cases and the PH distribution is also a best bet and approximation for a general distribution. Further the 
inventory systems with PH distributions are most general in forms and almost equivalent to inventory systems 
with general distributions. The bulk demand models because they have non-zero elements or blocks above the 
super diagonals in infinitesimal generators, they require for studies the decomposition methods with which 
queue length probabilities of the system are written in a recursive manner. Their applications are much limited 
compared to matrix geometric results. From the results obtained here, provided the maximum demand and 
supply sizes are not infinity, it is established that the most general model of the PH inventory system with bulk 
demand and bulk supply admits matrix geometric solution. Further studies with block circulant basic generator 
system may produce interesting and useful results in inventory theory and finite storage models like dam theory. 
It is also noticed here that once the maximum demand or supply size increases, the order of the rate matrix 
increases proportionally. However the matrix geometric structure is retained and rates of convergence is not 
much affected. Randomly varying environments causing changes in the sizes of the PH phases may produce 
further results if studied with suitable partition techniques. 
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